
STOP 



Early Journal Content on JSTOR, Free to Anyone in the World 

This article is one of nearly 500,000 scholarly works digitized and made freely available to everyone in 
the world byJSTOR. 

Known as the Early Journal Content, this set of works include research articles, news, letters, and other 
writings published in more than 200 of the oldest leading academic journals. The works date from the 
mid-seventeenth to the early twentieth centuries. 

We encourage people to read and share the Early Journal Content openly and to tell others that this 
resource exists. People may post this content online or redistribute in any way for non-commercial 
purposes. 

Read more about Early Journal Content at http://about.istor.org/participate-istor/individuals/early- 
journal-content . 



JSTOR is a digital library of academic journals, books, and primary source objects. JSTOR helps people 
discover, use, and build upon a wide range of content through a powerful research and teaching 
platform, and preserves this content for future generations. JSTOR is part of ITHAKA, a not-for-profit 
organization that also includes Ithaka S+R and Portico. For more information about JSTOR, please 
contact support@jstor.org. 



EXTENSION OF THE SIEVE OF ERATOSTHENES TO 
ARITHMETICAL PROGRESSIONS AND 
APPLICATIONS* 

Bv J. C. MOHKHEAD 

ALMOijiT without exception the methods einploved in the elaboration of 
the various existing factor-tables are modifications, more or less remote, of 
the classic Sieve of Eratosthenes. In most instances these factor-tables are 
designed to show the prime numbers and the prime factors of the composite 
numbers occurring in sets of consecutive integers. t The Sieve method, as 
applied to the determination of the primes in a set of consecutive integers 

m, m + 1, m + 2, m + 3, • • • m + n, 

may be briefly described as follows : We strike out from the numbers of the 
set the multiples of 2 ; from the remaining numbers of the set the multiples 
of 3 ; then, from those j'^et remaining, the multiples of 5, and so continue, 
for the multiples of all primes not greater than ^(m + n) . The numbers finally 
remaining in the set are all prime. If, as the procedure just described is 
carried out, the multiples of 2, 3, 5, • • • are conveniently indicated, we ob- 
tain what is called & factor-table for the set of numbers considered. 

The present pajjcr describes Sn analogous method applied to successive 
numbers in arithnjctical progression. The discussion in §§1-.t is quite 
general, while the remaining part of the paper is occupied with the discussion 
of applications and special cases. 

1. Extension of the Sieve method. Any arithmetical progression 
of m integral terms may bo written, in direct or inverse order, in either of 
the forms 

(1) «'■ + 6, 2«* + b, 3«* + b, ■ • ■ Hi«* + b, 

{!') — a* + b, — 2a* + b, — 3o* -^b, ■ ■ ■ — ma^ + 6, 

♦Head before tlie American Mathematical Society, April 29, 1905; December 29, 190S. 

t The Factor.Tables of J. Bnrckliardt, Z. Dase and J. Glaisher applying to all integers 
under 9000000, their recent extension by L. D. Lehmer, under the auspices of the Carnegie 
Institution, Washinaton, D. C, to 10000000, and the tables of "High Primes' of A. Cun- 
ningham and H. ,T. Woodall are the principal tables of the sort that have been satisfactorily 
verified. 

(88) 
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if suitable integral values are assigned to the symbols, a, b, k, m. The fol- 
lowing general discussion may be limited, therefore, to progressions of the 
forms (1) and (1'). Evidently vvc may assume, without loss of generality, 
thit a is positive and k positive or zero. We shall further assume, unless 
the contrary is expressly stated, that b is prime to a* ;* for each term in (1) 
and (1') would admit the common divisor of b and a'- as a factor, — a case 
which it is useless to include in a search for prime numbers or in the formation 
of a factor-table. 

The methods and results obtained in this paper have for basis and starting 
point the following 

Theorem. Let p be any prime number thit does not divide a'', and let the 
leant non-negative solution of the congruence 

xa^ + b = mod j>, 

for assigned values of a, b and k, be 

Then only the 

(^i-pY^, (P + ^L-pY^, {^P + ^kpy\ • • • ipm, + x,^y'' 
numbers in (1) are divisible by p and only the 

(P - xtpr\ (2i> - x^y\ {3p - x,,y\ . . . (pm'^ - x^y^ 

numbers in (1') are divisible by p, mp and m'p being the maximum integers 
such that 

pnip + Xkp, pm'p — Xtp ^ m.t 

As here defined, Xj.p is clearly a function of X", p, a and b. This will oc- 
casionally be indicated by writing x^pab- However, when it is believed that 
no ambiguity or confusion will arise the symbol will always be abbreviated to 
Xtp or even x^. 

The above theorem leads at once to the following 

• An exception is made in tlie application of tiie relation (il/) in §11. 

t From well known propt- rties of linear congruences it follows that no solntion xtp exists 
when p is a factor of a* ; bnt for all other jtrime values of p an unique solution exists such that 
^ Xt„ ^p — 1, and that all other solutions are represented by 

x = np+xi^, H=± 1, ± 2, ± 3, • • • . 

The Talldlty of the theorem is thus immediately rendered apparent. 
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ExTENSiox OF THE SiEVE OF Eratosthenes. Let V be the greatest prime 
not greater than ^M, where M is the maximum numerical value of the numbers 
in (1) and (1') for an assigned set of values of a, b, h and m. Let the values 
of x^p be known for all prime values of p us far as ir. Then to determine 
the primes in (1) we strike out from (1) the 

{x,^y\ (2 + x,.,y\ (2X2 + x^,y\ (3X2 + x^,y\ ■ . ■ (2m, + x^y^ 

numbers, which are multiples of 2 ; then the 

(^w)'", (3 + Xi.3)«", (2x3 + x^y'\ (3X3 + x,3)">, • • • (Smg + a:«,)"' 

numbers, irhich are luultiples of 3 ; fhe)i the 

(.r,,)'", (.-. + .r,,)'", (2 X :■) + •>■,.,)•". V> X a + x,,y\ ■ ■ ■ {bm, + x,,)"' 

numbers, which are divisible hi/ "» ; and so contiime, until fincdlij we strike 
out from (he remainimi numbers the 

(..•,,)•", {-jT + .T,„)'>', (27r + x,,„y^\ (;'>7r + .T,„)t", . . . (Trm, + x,.,)"' 

numbers, which are multiples of ir ; the numhers noio remaining are (he primes 
occurring in ( 1 ) . 

Similarli/, (o determine the primes in (1'), we strike out from (1') the 

(2 - .-i.,)'", (2 X 2 - X,.,)'", (3 X 2 - x,,,)t", • • • {2m, - Xi^^ 
(3 - X,,,)'", (^ X 3 - X,,)'", (3 X 3 - X,,,)"', . . . (Zm, - x,^)«" 

(T - x,,)t", (27r - x,„y\ {:W - x,,)'>' . . . (7r;«, - x^,)"' 

numbers, and (hose Ji nail i/ remuining in (1') are prime. 

If, as this process is carried out, the niuUijdes of 2, 3, 5, 7, • • • tt are 
indicated, we obtain a. fac(or-(able for (he progressions (1) and (1'). 

When Xj.,, x^.j, x,.}, • • • x,;, are known, this "sifting" is essentially iden- 
tical with the ordinary Sieve. In order, theretV)re, to render this extension 
of the Sieve method in any degree practicable, it is ne<!cssary to devise means 
of calculating expeditiously the numbers x,.^, for successive primes and for 
assigned values of a, b and A;. 

2. The calculation of x^p. Connection with residue tables. 

The choice of the most convenient means of computing x^^, i. e., of solving 
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the congruence 

(A) xa* + 6 = modp, 

depends, in general, on the relative values of k, a and p. The direct solution 
of (A), even for small p and for the few cases where residue tables are 
available giving the least remainders (mod p) of a*, is quite tedious. How- 
ever, when such tables are at hand, we express the solution of( ^) in the form 

(B) x,.p = — —J. = — ba"-'' = — fta"-*' mod p, 

where v is the principal exponent* of a (mod p), and k' is the least positive 
remainder of k (mod v), take the least remainder, R, of a''~^' from the tables, 
and compute directlj x^.^ as the least positive remainder (mod pi) of — bll. 
When i = ± 1, we may thus read the successive values of Xj.p directly from 
the residue tables. t Thus for a = 6, 6 = — 1, p = 115, k = 8, we have 
Xg = (>«-« = 7 mod 13. 

3. Calculation of x^-^i from Xf.. First method. In many cases, 
especially when series for diflPerent values of k, but fixed values of a and b, 
are under investigation, it is found expedient to employ one of the following 
methods of deriving ai^+Lp from Xj._p, x^.+j^^, from x^c+i.p^ "^"^ so on, be- 
ginning with a known Xf.j,. Setting A; = in (-6), we obtain 

(C) Xop=: — bmodp. 

Hence x^p is in every case known or readily determined, and therefore the cal- 
culation of Xj. for successive values of k may be begun with x,. 
Replacing x by x^, and writing (A) in the form 

^■a*+i + 6 = 0mod», 
a 

we have at once 

♦The rainimnm integer, ». snch that a" = 1 mod p. 

t The only printed residue table of suflicient extent to be of use in tlils connection is 
A Binary Canon, by Lt.-Coi. A. Cnnningbam. This tabic gives the least positive remainders of 
2' for x = 1, 2, 3, • • • »>, and for ali primes nnder 1000. An extension of tlie table (as yet in 
MS.), by Col. Connlngiiam and Mr. II. J. Woodall. includes all primes under 10000 for x g 100. 
and ail primes under 12000 for x ^ ."56. Hence when a =2, p < 12000, iind k is snch that v — ki 
falls within the range of values of x just stated, the Binary Canon is directly applicable to the 
computation of xi^. 
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where y is an arbitrary integer. Since a and j) are mutually prime, an unique 
choice of y from the numbers 0, 1, 2, • • • a — 1 is possible, such that 
(xj. + yp)la is integral and less than p. When y is so chosen,* we may 
write instead of (Z)) 

^k + yp 



(£) 



Xk+i = 



For example, when a = 2, y = or 1 according as Xt is even or odd ; 
and therefore the computation of x^.^-^ from x^ involves simply division b}' 2 
when x^ is even, or the addition of p and then division by 2 when Xt is odd. 
Thus, for a = 2, & = 1 and p = 37, Xo = 36 ; Xi = 36/2 z= 18 ; x„ = 9 ; Xj = 
(9 + 37)/2 = 23 ; x^ = (23 + 37)/2 = 30 ; etc. 

Again, when a = 10, y may be chosen at sight so that the final digit of 
yp + X/. will be zero.t Thus, for a = 10, i = 3, audp = 43, we have Xq = 40 ; 

..J = 4 ; x.> = (4 + 2 X 43)/10 = 9 ; x^, = (9 + 7 X 43)/10 = 31 ; etc. 

4. Calciilation of x^.,., from x^.. Second method. For brevity, 
set fj. = x^p„i and |j|. = Xt^„_,. The following relations are always satisfied, 
fi, ^j. and Xjt referring to the same a and j). 



(F) 

in 

(G) 
{G') 

(J) 
(J') 



X,. 



■i + l = — fl^i 
%• + ! = ^l^i 

^t + h = (-Ii)*a5i 

^k + h ^ ^l^k 
^k + h = ~ ih^k 
Xk + h = !»«* 
^k+nh— (— ^*)"** 
^k + nh = ^h'^k 



modjj. 



* The selection of y involves tlie solution of the congruence yp + x^=0 mod a, the diffi- 
culty of which, compared with tlie labor of solving (^4) direi-tly, will generally depend on the 
relative values of a, p and A;. Tlie application of {E, is essentially the same for A: small or 
large, which is not the case with the solution of {A). Also, if p > a, the application of {E) 
is clearly the preferable method. 

t in the case a ^ 10 the use of ordinary commercial addiny machines greatly facilitates 
the de ivation of a;^ , from x^; for the process reduces to the addition of p to x^ repeated until 
the Qiial digit of the siun is zero, and then dropping the cipher. This procedure may be modified 
so as to be applicable to certain other special cases. 
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Evidently (F), {G) and (H) are contained in (J), and (F'), (G') and 
(i/') are contained in (J'). All, however, are useful in computation, and 
are written together for convenient reference. 

The demonstration of (J) and (J') establishes the validity of the six 
pi'eceding relations as well. From the definitions of x^. and 1^ we have 

(c) X,. = — ba~'' mod p, 

(7 ) — f ft = a-* mod p. 

Multiplying ('•) by (7) n times, member for member, we obtain 

(_ f^)»xj = - Ija-<-^+>^'> mod p. 



But by definition 
therefore 



a-j ^. „ft = - 6«-(* + "*' mod p ; 
%• + »* = (- ^h)"^-kT^odp. 



Tlie demonstration of («/') is similar. 

The computation of x^+i, Xt+ft, or x,.^„h from Xg. by means of the above 
relations involves only the direct calculation of the least positive remainder 
(mod J)) of a product. Thus, forp = 41, o = 14 and b = 5, we have, from 
(C) and (F), 

^i=l, |( = (1 + 41)/14 = 3; Xo = 36, 

and therefore, by (F'), 

Xi = 3 X 36 = 26 
a;. = 3 X 26 = 37 
X3 = 3 X 37 = 29 
x^ = 3 X 29 = .5 



mod 41. 



Thus (F) and (F') are admirably adapted to the compution of x^ for succes- 
sive values of k, yielding a method generally practicable. For a few special 
values of a, however, the method of §3 involves less reckoning and is there^ 
fore preferable. Any one of the formulas (G), {G'), • • • (J') may be em- 
ployed in calculating Xf. for non-consecutive values of k, — in practice, that 
one is selected with which the smallest | or ^' may be used. Applied in this 



94 HOHEHEAO 



[January 



way at intervals, the congruences (G), (G'), • • • {J') are especially valua- 
ble as check formulas. 

5. Relation between x^-paj, and x^pa^b- For brevity, we write 
^k = ^kpab and x'f. = a;ipa_6. By definition 

aji-a* + b = mod p, 

a;[.a* — ft = mod p. 
By addition, we obtain 

(Tj, + xj.) a* = mod^ ; 
hence 

a",, + aij. = mod ^, 

and since x^, a;^ < j), we have at once 

(K) Xk + x'^=p, 

except when Xf. = x',. = 0, i. c, when b is divisible by p. This exceptional 
case, 6 divisible by p, requires no computation, and may therefore be excluded 
from discussion in this and the following section.* 
If b is replaced by — &, formula ( C) becomes 

{C) x'op = b mod jj, 

which could have been inferred also from (/i") . It follows as an immediate 
consequence of the definition of x'^., that the same methods and rules that 
apply to the computation of x^. apply equally well to the computation of x[., 
starting from the initial value given by ( C") . 

When one of the numbers x^., xj, is known, the other may be obtained by 
the relation (K ) ; however, when a table of successive values of Xj. and x^ 
together is being formed, as illustrated in §6, (K) is employed most effect- 
ively as a check formula. 

6. Calculation of Xj. and x[. when a = 2. As has been noted in 
§4, the calculation of x,,+ , from .%. by means of formula (E) is extremelj- 
simple when a = 2, involving onl}' division by 2 when Xj. is even, and the 

• The relation <K) is contained implicitly in the theorem of Jl, as may easily be seen if 
we note that the numbers in the series (1) and (l')> with the signs changed throagbont, refer 
to — 6. In fact, the x« for the numbers in (1'), with the opposite signs, is simply p — «*, t.f., i'». 
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addition of p and division by 2 when x^. is odd. Even greater simplicity is 
attained when the computation of successive values of x^ and xi. is performed 
simultaneously by a method based on the relations :* 

(X) ^k -i- ^'k + i = ^* + i' ^''i^° ^k is odd, 

(L') x'^ + %.+ , = xit^.], when xj. is odd. 

We shall demonstrate only (L), as the demonstrations of (L) and (L') are 
similar in form. 

By definition of Xj, we have 

x^a* + 6 = mod^, 
and multiplying by 2, 

(k) Xj,«* + i + 26 = modjp. 

Also, by definition of xi + i, 

(k') xJt + 1 a* + > - 6 = mod j9. 

Adding (k) and (k') , we obtain 

(.Tj + xj.^i)«*' + i + 6 = mod 7?, 
and therefore 

(?) xt + x[. + i = Xi. + imodp. 

In keeping with the conventions stated in §1, ^> = 2 is excluded when a = 2, 
and hence it follows from (/i ) that one of the numbers x^, xj. is even and tiie 
other odd. Let x^. be odd ; then xj. is even and xj.^.i = (Xi.)/2. But, from 
(K), Xi + x'l, =p, and hence x^. + 4 + i = x^. + (xi.)/2 <;). We may there- 
fore, when Xfc is odd, replace the congruence (/) by the equation 

The formulas (L), (L') and (E) yield the following compact method for 
the calculation of successive values of x,. and xjr. when a = 2. 



* (i) and (V) are valid only for a = 2. The generaliisatlons take the forms 

(m) Xj + (a - l)r[ ^ , = z^ ^ , mod;), 

Cm') a;^ + (« - I)x, ^ , = xj ^ J m..dp, 

respectively. These congniences are not especially nsefnl in computation wlien a > 2, since 
when one of the numbers a;^ + f ^J + 1 1» Isnovyn, the other can be found from {K ) more readily 
than from (»») or (m'). 
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The numbers aio and x'^ are known and one is even ; hence one of the 
numbers Xi, x[ can be obtained by a single direct division by 2, and the other 
bv the application of (L) or (L'), involving the addition of two numbers, 
each of which is less than p ; then one of the numbers x^, a;^ can be obtained 
by a single direct division by 2, and the other by the application of (i) or 
{L') ; and so on, for X3, 3-3; X4, x'^; • • •. 

The compactness of the method and the speed it is possible to attain in 
the reckoning are perhaps best understood in connection with a few examples 
arranged in tabular form as follows. The number heading each column rep- 
resents k, and the upper and lower numbers of the pairs in the columns rep- 
resent the corresponding values of x^ and x^ respectively. 



k = 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


p = 67|.r,= 
h= 7j4 = 


60 

7 


30 
37. 


15 
52 


41 
26 


54 
13 


27 
40 


47 

20 


57 
10 


62 

5 


31 
36 


49 

IS 


p= 16231 1.rj = 
6= 127)x^== 


16101 
127 


8052 
8179 


4026 
12205 


2013 
14218 


9122 
7109 


4561 ll039S 
11670 5S35 


5198 
11033 


2599 
13632 


9415 
6816 


12823 
3408 


p = 431 |.r, = 


430 

1 


215 
216 


323 
108 


377 
54 


404 
27 


202 101 
229 j 330 


266 
165 


133 
298 


282 
149 


141 
290 



Li the above table, the numbers in the lows Xj. or the rows x',. that are 
obtained by division by 2 are sufficiently evident. As illustrations of the 
application of ( L) and {L'), we have, for^ = 67 and 6=7, Xj = 7 4- 30 = 37 ;. 
xl, = 37 + 1.5 = ,52 : .T3 = 1.5 + 20 = 41 ; cc^ = 41 + 13 = 54 ; etc. For^^ = 
16231 and b = 127, .r„ = 4,5(;i + 5835 = 10396 ; x'- = 5835 + 5198 = 11033 ; 
etc. 

As illustrated in the above cases, we alwaj's have, when « = 2 and b is 
numerically less than p, Xj = (p — b)/2 and x{ = (p + b)/2, 

In the manner indicated one may, after a little practice, compute a table 
of the kind indicated above almost as i*apidly as the numbers it contains can 
be written down. For the single division by 2, or else the single addition, 
by which each of the numbers that are to be entered in the table are obtained, 
may be performed mentally as the numbers are recorded. Furthermore, when 
Xj. and x[. are calculated by this means, the formula (fT) yields an independent 
and effective check, and one that can be applied at a glance. When the ar- 
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rangenieat of the table is as above, for example, the sum of the two numbers 
ia each rectangle gives the corresponding prime. However, it is unnecessary 
to check each pair, x^., x-j., but it is sufficient to apply the test at intervals of 
a half score terms or more. If ordinary care is exercised, a table thus com- 
puted and checked needs no further verification. 

7. A table of values of x^ and x[. for « = 2.* The author has 
undertaken the elaboration of a table of values of x^. and xj. for a = '2,b = 1, 
^• = 1, 2, 3, . . . 100, and for all primes under 100000. This table is now 
under way. A table of Xj. alone for A; = 1, 2, • • • 19 and for all primes un- 
der 25000 was computedf by the method of §3, but the formation of the 
larger table for both Xj. and x[. was subsequently begun anew by the more rapid 
method of §(>. The arrangement of the table is that shown in §6, except that 
the first column is headed ^j and contains in order the successive primes 3, 5, 
7, 11, • • • and the column for k — is omitted. 

The table, when completed, will be applicable to the determination of 
primes and the formation of factor-tables of numbers of the forms >i2'^ i 1 up 
to 10"*, and to the determination of all factors under 100000 of numbers of 
these forms beyond 10'", for^• S 100. 

8. An application of the table. We pass now to the application 
of the table to the determination of the primes in the two arithmetical series : 

(2) 219 ^ 1^ 2(2i») + 1, 3(2i») + 1, . . • 1211(219) + i^ 

(2') _2»9+l, - 2(219) + 1, _3(2'9) + 1, ..._ 1211(219) + 1. 

If we include 1 with the numbers in (2) and (2'), we may write instead of the 
two seiies the single arithmetical series : 

(2") - 1211(219) + 1, - 1210(2'9) + 1, • • • 

-219+1, 1, 219 + 1, 2(2'9) + 1, ... 1211(219) + 1. 
Writing down in order the coefficients of 2'9, we obtain 
(II) - 1211, - 1210, -1209, . • • 

- 3, - 2, - 1, 0, 1, 2, 3, . • . 1209, 1210, 1211. 

* A paper of Mr. f^. L. Dines, to be published In the April number of the A.nnals, con- 
tains a detailed discussion of the case a = 6 and & = 1, with various applications and the de- 
scription of a table he has elaborated. 

t By the .author and an assistant, Miss C. Ame«, in April, 1905. This table is Included in 
the appendix of the author's thesis on •' Fermat's numbers and numbers of the forms 2*3 ± 1" 
on file in the Yale University Library. 
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Referring to the column headed lit in the table, we find 

•"^ni.H = 1> •''19,5 = "» •*'li),T == "'' •'-19,1] =.',••• 

Striking out from the numbers in (II), the positive and negative multiples 
of 3, increased by 1, —namely, 1, 4, 7, 10, 13, • • • 1201, 1204, 1207, 1210; 
-2,-5,-8,-11,-14, . '■ . - 1202, - 120o, - 1208, - 1211 ; then from 
the numbers remaining in (II) the multiples of 5, increased by 3, — namel}', 
3, 8, 18, • • • 1203, 1208 ; - 7, - 22, ... - 1207 ; then from those remain- 
ing the multiples of 7, increased by 3, — namely, 17, 24, • . • 1200; —4, 
— 18, ... — 1194 ; and so continuing, for all primes under \/\ 1211(2^*)+ 1 (, 
the numbers finally remaining in (II), except 0, are the coefficients of the 
prime numbers that occur in tlie series (2"), or, which amounts to the same, 
in (2) and (2'). 

The author has cari-ied out the work for this series as indicated, and has 
found that the following numbers in (II), and these only, correspond to 
primes in (2") : 

1198, — 119S, —1186, — llf>4, —1156, — 1U3, —1141, —1134, —1129, —1126, —1125, —1123 
1108, — 1104, — 1095, — 1090, — 1084, — 1083, — 1060, — 1045, — 1044, — 1039, — 1030, — 1021 



-978, 


— 973, 


— 958, 


±936, 


— 925, 


— 924, 


— 919, 


— 910, 


± 8!»4, 


— 888, 


— 885, 


— 883 


- 808, 


— 853, 


— 850, 


— 828, 


— 804, 


— 801, 


— 793, 


— 778, 


— 759, 


— 754, 


— 748, 


— 745 


-735, 


— 733, 


— 730, 


— 720, 


— 705, 


— 700, 


— 696, 


— 685, 


— 664, 


— 663, 


— 654, 


± 609 


- (103, 


— 001, 


— 600, 


— 593, 


-579, 


— 550, 


— 540, 


— 534,, 


— 520, 


— 493, 


— 484, 


— 463 


- 460, 


— 450, 


— 444, 


— 441, 


— 421, 


— 415, 


— 414, 


±411, 


— 408, 


— 400, 


— 394, 


— 391 


b 369, 


— 363, 


— 360, 


— 339, 


± 324, 


— 315, 


— 313, 


— 306, 


— 303, 


— 286, 


— 268, 


— 265 


b264. 


— 241, 


— 234, 


-229, 


— 223, 


— 213, 


— 208, 


— 199, 


— 153, 


— 150, 


— 139, 


— 133 


-115, 


-114, 


— 103, 


— 99. 


— 93, 


— 91, 


— 70, 


— 75, 


-73, 


— 70, 


± 69, 


— 61 


— 51, 


— 36, 


-.34, 


— 1, 


















11, 


14, 


26, 


27, 


41, 


44, 


50, 


51, 


57, 


± 09, 


71, 


77 


131, 


132, 


134, 


147, 


150, 


I.VJ, 


179, 


191, 


194, 


197, 


200, 


212 


216, 


221, 


225, 


230, 


2,")2, 


200, 


± 204, 


270, 


282, 


290, 


302, 


309 


812, 


320, 


±324, 


335, 


341, 


350, 


365, 


± 369, 


380, 


390, 


404, 


±411 


417, 


419, 


422, 


435, 


440, 


447, 


4.-)9, 


467, 


470, 


476, 


485, 


489 


492, 


497, 


61(i, 


524, 


551, 


554, 


55.-), 


582, 


585, 


594, 


±609, 


627 


629, 


630, 


641, 


0112, 


671, 


680, 


090, 


714, 


719, 


720, 


732, 


787 


740, 


747, 


755, 


762, 


771, 


782, 


792, 


795, 


809, 


810, 


819, 


827 


837, 


854, 


875, 


876, 


879, 


884, 


± 894, 


896, 


929, 


±936, 


942, 


945 


951, 


987, 


989, 


1001, 


1014, 


1019, 


1031, 


1049, 


1059, 


1077, 


1089, 


1094 


1100, 


1110, 


1116, 


1127, 


1136, 


1140, 


1154, 


1155, 


1161, 


1167, 


1169, 


1181 


1185. 

























The numbers above that are marked with a double sign, as ± 93C, occur 
in both the negative and positive coefficients. Whenever such is the case, 
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the corresponding numbers in (2"), as ± 936(2'®) + 1, differ numerically by 
2, and compose what is called a " prime pair." There are eight such pairs 
occurring in (2"). 

After the table had been made, the time occupied in applying it to the 
determination of the primes occurring in (2) and (2') was about eight 
hours. This is a shorter period than would be required in multiplying out 
the 2422 numbers in (2) and (2') and consulting factor-tables like those of 
Burckhardt, Dase and Glaisher, if such factor-tables were extended as far as 
634912 769, the largest number considered. All the numbers in (2) and (2') 
in which the coefficient of 2'* is numerically greater than 19 lie beyond 10^ 
and hence beyond the range of the factor-tables. 

In the above case, while the numbers in (2") range in numerical value 
from 524 287 to 634912 769, the "sifting" process is applied simply to the 
coefficients of 2^*, a set of consecutive integers i"anging from — 1211 to +1211. 
The power of the general method of this paper is due to just this feature. 
After a table of x^.p of sufficient extent to meet the reciuirements of the case 
in hand has been computed, the sifting is applied not directly to the numbers 
under investigation, but to a range of smaller consecutive integers, the coeffi- 
cients of a*. Thus is obtained a method, not difficult to apply, for the deter- 
mination of the primes and the factorizatiim of the composite numbers of a 
given form xaJ' + b, where the investigation of each number separately would 
be quite impracticable. 

9. Application to the factors of Fermat's Nxunbers. As is 
well known, all factors of Fermat's Numbers,* i^„ = 2*" + 1, must be of the 
form 

(?„ = y.2'' + '^+ 1. 

All the known factors of F„ were identified by testing the divisibility of /'„ 
by primes of the form Q„. The primes of the form 5'. 2'* + 1, given in §8, 
were tested as possible Fermat factorsf for q = 11, 14, 26, • • • 989, with the 
result that no Fermat factors were found in addition to those already known. 
Hence, for u > 16, there are no more factors of F^ in addition to those pre- 
viously known, under 1001. 2i» + 1, or 524812289. 

• Fn is known to be prime for n = 0, 1, 2, 3. 4, and composite for n = 5, (>, 7, 9, U. 12, 
18,23. 36, 38, 73. See Proceedings of the London Mathematical Sorietii, ser. 2, vol. I, p. 175; 
vol. 3, p. XXI; vol. 5, p. 237; Messenger of Mathematics, vol. 37, p. 65; Bulletin of the American 
Mathematical Society, vol. 11, p. 343; vol. 12, p. 449. 

t 7. e., factors of Fermat's numbers. 
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Bj making use of the 55"" and To**" columns of the table described in §7, 
in its present extent, it was found that neither 5 . 2^ + 1 nor 5 . 2'^ + 1 admit 
factors under 2000. Subsequent tests showed that both are prime, and that 
5 . 2" + 1 is a factor of F-y The last result was obtained hy the following 
method. 
Setting ^ = a . 2'* + 1 , we have 

5 . 2'3 + 1 = modj:>, 

and hence, using the notation of §4, 

664 = ^Mp21 = 5.2. 
Then applying (//) we obtain successivelv 

I,, = - ^^ = _ 5-^ . 2» moAp 

fj, = -b* . 2« 

^3. = - i)« . 2«« = 5' . 21- 

f.,« = 52« . 10^. 

From this stage on, the numerically least remainder (mod^) of each ^\n was 
found by direct division by ^;. This process showed finally f^-, = 1, and there- 
fore that 2-" +• 1 is divisible by p. Thus the identification of the Ferniat 
factor 5.2" + I was eft'jcted solely by the methods of the present paper,* and 
less arduously, as far as f^.u, than by the usual method of calculating succes- 
sively the remainders (mod jt) of 2'-', 2*", • . • 2-"'. 

A more detailed report on the application of the methods of this paper 
to Fermat's factors will be held, for the sake of completeness, until the table 
of §7 is further advanced. 

10. Connection with Lebon's table. Professor Lebon, in a paper 
published in the Bullelin of the Amei'kan Mathenintical Socktif, vol. 13, 
p. 74, describes a process for the construction of a fable ai)plicabie fo the fac- 
torization of the numbers in arithmetical progressions in which the common 
difference is the simple product of consecutive primes, 2. 3. 5. • • • p^, and 
the first term is any one of the <^(2. 3. 5. • • • p^) numbers less than and 

• Tills method is stated briefly in dilfereut form iu Bulletin of the American Mathemati- 
cal /Society, vol. 12, p. 237. 
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prime to the common difference. The common difference and first term, 
called the base and index of the progression, are represented by £ and / rr- 
spectiveh'. The system oi (f){B) arithmetical progressions whese general 
term is BK + I, where K is successively equal to the positive integers starting 
from zero, and / is prime to B, contains all prime numbers except 
2, 3, 5, . • • p^, and all composite numbers except multiples of 2, 3, 5, • • • p^. 
The numbers recorded in the table are the mininmm values of Ii that 
render BK + /divisible by each successive prime following /);^.* Evidently 
a table so constructed is applicable to the decomposition of, or the deter- 
mination of the primality of successive numbers in progressions of this special 
form, as far as pi (the last prime included in the table), or to the inves- 
tigation of the divisibilitj' of isolated numbers in such progressions, in a 
manner similar to that described in the present paper. In fact, if we take 
a = B, 6 = 1 and h = 1, the numbers »i occurring in my table in the first 
column following the primes coincide, for the primes following p^, with the 
minimum values of K recorded in Professor Lebon's table. If the latter 
table is extended to include the <i>{B) values of 7, the complete table is 
clearly applicable to successive integers, aside from multiples of 2, 3, • • • jjj, 
and hence applicable to the formation of factor-tables of the usual form. 

However, the method of Professor Lebon is not so easily and directly 
applicable as la the method of the present paper, to the investigation of the 
divisibility of numbers of the form w««* + b when k exceeds unity. Thus the 
special form of the numbers in. 2'- + 1, treated in §§«, 9, would be of no ad- 
vantage if the table of Professor Lebon were to be applied. 

11. Application to ranges of consecutive integers. The method 
of this paper may be extended so as to be applicable to a range of consecutive 
integers without restricting the form of the progressions as in §10. Let a 
and k be fixed. Then, like the Lebon progressions, the system of 4>(a^') 
arithmetical progressions represented by m«^ + 6, where m = 0, 1, 2, • • • and 
b assumes successively the <^(a*") positive integral values less than and prime 
to a*, contains all positive integers that are prime to a. Consequently, a table 
of values of x^.p^ for the primes as far as tr and for the^(a*) values of b would 
suiBce, like the Lebon table, for the decomposition, or the determination of 
the primality of all numbers prime to a as far as tt^. 

* A table compiled by Professor Lebon lias been published under the title Table de 
caractiristiques relatives a la base 2310 des faeteurs premiers d'lin nombre inferior a 30030, 
Paris, 1906. 
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The efficiency of such a table, as well as the amount of labor involved 
in the compilation, will be determined by our choice of a and k. To fix ideas, 
suppose it is desired to construct a table eflfective to 10*. A convenient 
choice would be « = 10, k = 4. Then the system of ^(10*), or 4000, pro- 
gressions of 10000 terms each, whose general term is m.lO* + b, contains 
every odd number except multiples of 5, from 1 to 10* — 1. To be effective 
as far as 10*", the corresponding table of values of Xtp^g, should be calculated 
for 1226 primes, viz., 7, 11, 13, • • • 9973.* The following is a convenient 
arrangement of the table. Here a;4pio6 stands at the intersection of the row^ 
with the column b. 



a = 10, k = 4. 



h = 


1 


3 


7 


9 


• 


p = l 


.5 


1 





3 


• 


11 


10 


8 


4 


2 


• 















233 



237 



997 



432 



299 



33 



897 



956 



690 



To illustrate the method of applying the table, suppose we wi,sh to 
decompose the nninbor 26840 237, = iV, say. It is easil}' seen that 3 is not a 
factor of JV. Writing X= 2684 X 10* + 237, and referring to the column- 
headed 237, we find 

26t<4 # 2 mod 7 

5t 5 mod 11 



= 690 mod 997. 



Here 690 is the ./?)•.«« remainder obtained on dividing 2684 bj' a prime that 
coincides with the corresponding x^j. Hence 997 is the least factor of j\L 
The remaining factor, 26921, is less than 997* and therefore prime. Thus we 
have the complete decomposition, iV= 997 X 26921. If for the primes as far 



* It is useless to include 3, as divisibility by 3 is easily tested. 
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as y/iVno remainder were found to coincide with the corresponding ajpj,, then 
would iV be prime. 

The application of the table would be greatly facilitated if a residue table 
were available, giving the least positive remainders of each integer as far as 
10* for the prime moduli 7, 11, • • • 9973. If the arrangement of the residue 
table were similar to that of the above table, then to effect the decomposition 
of the number 26840237, for example, we should need only to compare cor- 
responding numbers in the columns headed 2684 and 237 in the residue and 
factorization tables respectively, until a coincidence is discovered at^ = 997. 

The application of the above table has one evident advantage over that of 
Professor Lebon, covering the same range of numbers. We recognize at 
sight the coefficients required to express a number in the form m.lO* + 6; 
while a division must be performed to express a number in the form 
m('2 . 3. • • • p^) 4- b preliminary to applying the Lebon table. Aside from 
this feature, the two tables differ little in efficiency when a residue table is 
available with each. 

The labor of compilation is not so great as to render the task inipi-actic- 
able. Starting from the value of x^ given by (O), §3, we calculate succes- 
sively, for given values of h and p, Xj, x.^, x-^ and finally 054 by the method of 
§3. Aided by the use of calculating machines, a computer of average speed 
can evaluate by this method and record about one hundred and twenty x^'s 
y>er hour, when the corresponding primes contain not more than four places. 
The labor may be further reduced by means of a relation which we shall now 
derive. 

Let us write, for brevity, x,, for x^|,l^|,. We have, by definition, 

10* ■ X,, + b = modj) 

and 

10* . X,,. + // = mod jj. 

Wiicnce 

10*(x,, 4- X,,.) + f/ + b' = niod;j, 
and therefore 

(^f) XV, + ;.' = X,, + X,,. mod p. 

When the numbers in colunni 1 have been found by the method of §3, 
repeated applications of (31) will give all the remaining coimnns. However, 
the relation (if) can l)e used most efFectiveh' if two auxiliary columns have 
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been calculated, viz., for 6 = 2, 4.* Since two successive values of b in 
the table differ by 2 or 4, the addition of x^ or x^ to any x^ and the sub- 
traction of p if the sum exceeds p — 1 will give the corresponding a;^ .(- ^ or Xj, ^. 4 
in the next column. Thus for p = 997, we have x^ = 864, x^ = 731 ; 
Xj, = a;, + x.2 = 33 + 864 = 897 ; Xjjg, = x^gg + x^ = 956 + 731 = 690 mod p. 

If a calculating machine is used, a computer can evaluate by means o{{M) 
and record from four to five hundred X(,'s per hour, for primes of four places 
or less. The table may be checked by applying the method of §3 to every 
tenth or twelfth column. f We may estimate, therefore, that a single com- 
puter would require four to five years for calculating, recording and checking 
the Xj's for the entire table of 4000 columns and 1226 primes.} 

The generalization of the relations developed in this paper for composite 
moduli, and the incidental application of the methods of deriving aj,p„|, and 
^'ipah to the calculation of power residues are only remotely connected with 
the purpose of this paper and will be reserved for later publication. 

NOHTHWESTBRN UNIVEUSITY, 
EVASSTON, IlX. 

APKIL, 1908. 

♦ Here h is not prime to a, so that these two columns should not be Included in the table. 
It is easily seen that (M) is valid for ft = 2, 4. 

t It is evident from the definitions that Xi = % ^ ,,. Hence the nnmbers in any row of 
the table recur in exact order at intervals of lOp in the values of 6. For primes under 1000 
this property facilitates the calculation to some extent. 

J The manuscript factor-tables of Professor Kulik, left In charge of the Vienna Academy, 
extend to 100330201, — a little beyond the proposed limit of the factorization table here de- 
scribed. However, the Kulilc tables are at present almost useless on account of the many errors- 
they contain. See Bulletin of the American Mathematical Society, vol. 14, p. 10*!. 



